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Pion form factors in two-flavor QCD S. Hashimoto"' 

1. Introduction 

Electromagnetic and scalar form factors of pion are fundamental quantities in the low energy 
dynamics of pions. In the chiral perturbation theory (ChPT) their radii are related to the low energy 
constants at 0{p'^). The electromagnetic charge radius is experimentally measured rather precisely, 
{r^)y = 0.452(10) fm^ [jl]], and the scalar radius can be related to Fk/Fji or nn scattering ampli- 
tudes. Therefore, these simple quantities provide a good testing ground of the lattice calculation 
techniques. In particular, the chiral extrapolation seems to induce large systematic errors unless 
the unquenched lattice simulation can treat pions well below 300 MeV where the chiral logarithm 
effect should become prominent [Q, ||]. It is interesting to check if the lattice data reproduce the 
expected chiral logarithms in the pion form factors. 

Our numerical calculation is performed on the two-flavor gauge ensembles produced by the 
JLQCD collaboration iQ] at j3 = 5.2 on a 20^ x 48 lattice with the 0(a)-improved Wilson fermions. 
High statistics is the key for the form factor measurements; we accumulated 1,200 gauge configu- 
rations for each five sea quark masses. The pion form factors are obtained as a by-product of our 
calculation of the Ki^, form factors [^]. It turned out that the largest effect on the Kij, form factor at 
zero momentum transfer ff^{0) comes from the shift of the zero recoil point to t,„ax = ("^i? — "^;r)^- 
For the evaluation of f+^{0)/f^^{tmax) we need a precise knowledge of the pion/kaon charge ra- 
dius, which is a subject of this paper. 

2. Pion form factors 

The pion electromagnetic form factor Gji{q^) is defined as 

W)\JMp))=G,{q^){p+pr. (2.1) 

where q^ = {p — p'Y and Jem is the electromagnetic current. Because of the charge conservation, 
the form factor Gji{q^) is normalized as Gji{0) = 1. It is well known that the experimental data 
in the small momentum transfer = —q^ < 1 GeV^ are well described by the vector meson 
dominance hypothesis Gji{Q^) = 1/(1 + Q^/m^). (For a summary of experimental results, see for 
example ||^.) 

The charge radius (r^)y is a slope of Gjt{q^) near the zero momentum transfer 

G^{q^) = l + ^{r^)^q^ + ---. (2.2) 
Near the chiral limit, ChPT predicts the chiral logarithm 
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(2.3) 



Unlike the pion decay constant, for which the chiral logarithm appears in the form m\\nm\, the 
chiral limit of the charge radius is divergent, and one expects a large effect as the physical pion 
mass is approached. The current experimental value is {r^)y = 0.452(10) fm^ [|l]]. 
The scalar form factor Gs{q^) is defined in a similar manner. 

{K{p')\qq\K{p))=Gs{q^)^ (2.4) 
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Figure 1: Double ratio (3.1) for the vector form factor at K = 0.1340 (left) and 0.1355 (right). Momentum 
is (1,0,0) (circles), (1,1,0) (squares), and (1,1,1) (diamonds) in units of 2n/L. 



Here the scalar operator is understood as a flavor non-singlet one. The scalar radius (r^) J is defined 



for Gs{q^) like (r^)y for Gji{q^) in (2.3). The chiral logarithm is even stronger for the scalar 
radius than for the charge radius, i.e. the coefficient of the log term is — 15/2(4 tt/)^ rather than 
-l/(47r/)2in(0). 



Although there is no direct physical process induced by (2.4), the scalar radius can be extracted 
from the tttt scattering using ChPT, (r^)J = 0.61(4) fm^ [Q]. It can also be related to /k/ fn, 
from which one obtains 0.65(4) fm^. For non-degenerate quark masses, the scalar form factor 
corresponds to f§^{q^) of the Ki^ decay. Its radius (r^)f'^ is measured as 0.235(14)(8) fm^ or 
0. 165( 1 6) fm^ [|l^] for charged and neutral kaons, respectively. It indicates that the flavor symmetry 
breaking is a large effect as suggested by the large chiral logarithm. 



3. Lattice calculation 

For the calculation of pion form factors we calculate three-point functions C^^''^\p7c,tj) with 
operators J (= V4 or S) inserted at time tj. Smeared pion interpolating fields are set at time t = 
and r/2 {T is the temporal extent of our lattice T = 48). Spatial momentum is inserted att = T /2 
and tj, so that the initial pion propagating from f = to is at rest and the final pion between tj 
and r/2 has a momentum p^^. To extract the form factor we use a double ratio 



C 



{P7t,tj) 



CP'JP\0,tj) Gj{q^)m„+EM 



C'''{Pn,tj) 



Gj{0) Irrin 



o<^tj<. r/2, 



(3.1) 



where the two-point function C^'^^^ {p^,t) in the denominator is constructed with the same smeared 
operator and a local operator P^. In the large enough time separation the numerator gives the 
ratio of form factors Gj{q^)/Gj{0) up to some kinematical factor, while the denominator becomes 
identity because (7r(/7)|P^|0)/(7r(0)|P^|0) = 1 must be satisfied for Lorentz invariance. Positive 
and negative tjs are averaged in order to increase statistics. 
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Plots are shown in Figure [T] for the vector current / = V4. Results for the heaviest (K = 0. 1340) 
and lightest (K = 0.1355) quark masses are plotted; the data have similar quality for other quark 
masses. (Sea and valence quark masses are set equal in our calculation setup.) We find a fairly 
good plateau beyond f ~ 4, though the signal for the largest momentum (1,1,1) is not satisfactory. 
We fit the data with a constant in the range [4,8] to extract the form factor ratio in the RHS of (3. 1). 



4. Form factor results 

In Figure ^ we plot the vector form factor for each K value. We fit the data with two different 
fit forms: 

1 



G7t{t) = - ^-^ + Cl(r2^)^ (free pole), (4.1) 

1 



Gjiit) = -^ + do{rit), (measured pole). (4.2) 

1 —t/my 

Here, t = is normalized by the Sommer scale tq. cq, c\, and do are fit parameters. Both forms 
are motivated by the vector meson dominance ansatz, but the pole mass is a free parameter in the 
"free pole" form ( [4.1[ ), while the vector meson mass measured on the lattice at a given quark mass 
is used in the "measured pole" form (42). In Figure ^the fit curves with the "free pole" form are 



plotted, but both can fit the lattice data equally well. 

Theoretically, we prefer the "measured pole" form, because the vector meson pole must appear 
with the known vector meson mass from the analyticity. There are other contributions from higher 
resonances and continuum states, that we may parametrize by the linear term do{rQt) near t = 0. 
This statement applies when the valence quark mass is heavy enough that the p —>■ nn threshold 
does not open, which is the case in our lattice calculation. As the quark mass decreases, the nn 
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Figure 3: Chiral extrapolation of the pion charge radius. Results from the "free pole" (ckcles) and "mea- 



sured pole" (diamonds) are plotted together with the extrapolation curves with ( 14.3| ). An expectation from 
the pole dominance model 1/my is also plotted (triangles). 



continuum state starts to contribute and finally it reproduces the ChPT prediction near the chiral 
limit. Therefore, in the physical quark mass regime the pole dominance form should be modified 
appropriately (see, for example, [pl|]). 

Chiral extrapolation of the pion charge radius is plotted in Figure ^. The "free pole" fit results 
have essentially no dependence on m| within the large statistical errors. With the "measured pole" 
fit the statistical error is much reduced and we find an upward trend toward the chiral limit. It is 
strongly correlated with the pure pole dominance model, which gives 6 /my and the upward trend 
is obvious. 

For the chiral extrapolation we employ the one-loop ChPT formula with a higher order analytic 

term 

1 

^'-'^v = Co - In +Cim2, (4.3) 

where Co and C\ are fit parameters and \l is an arbitrary scale. We can see that the chiral logarithm 
enhances the charge radius near the chiral limit, and the value extrapolated to the physical point is 
(r^)^ = 0.396(10) fm^, which is about 12% below the experimental value. 

For the scalar form factor we do not expect the dominance of the lowest-lying resonance (the 
ao meson), as it is already rather far from the = point. We therefore simply fit the form factor 
by a polynomial. Chiral extrapolation of the scalar radius is shown in Figure ^, where we draw 
the fit curve with the chiral log form plus a C?(m|) analytic term. The chiral logarithm is very 
strong for this quantity and the extrapolation could be sensitive to the details of the fit function. 
Our preliminary result neglecting such systematic effect is 0.60(15) fm^. 

This work is supported by the Large Scale Simulation Program No. 132 (FY2005) of High 
Energy Accelerator Research Organization (KEK), and also in part by the Grant-in-Aid of the 
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Figure 4: Chiral extrapolation of the scalar radius. 
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